Chandra Observations of Nuclear Outflows in the Elliptical Galaxy 

NGC 4552 in the Virgo Cluster 
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ABSTRACT 

We use a 54.4 ks Chandra observation to study nuclear outflow activity in NGC 4552 (M89), 
an elliptical galaxy in the Virgo Cluster. Chandra images in the 0.5 — 2keV band show two 
ring-like features ^ 1.7 kpc in diameter in the core of NGC 4552, as reported previously by 
Filho et al. We use spherically symmetric point explosion shock models to argue that the shape 
of the surface brightness profile across the rims of the rings and the temperature of hot gas 
in the rings are consistent with a Mach 1.7 shock carrying mean mechanical power Lghock ~ 
3 X 10'*^ ergss"^ produced by a ^^ 1.4 x 10^^ ergs nuclear outburst ~ 1 — 2 Myr ago. We find the 
gas temperature in the central ^ 100 pc of the galaxy to be 1.0 ± 0.2 keV, hotter than elsewhere 
in the galaxy, suggesting that wc may be observing directly the reheating of the galaxy ISM by 
the outburst. 



Subject headings: galaxies: clusters: general 
galactic medium - X-rays: galaxies 

1. INTRODUCTION 

X-ray observations of nearby galaxies provide 
a unique opportunity to study directly the inter- 
action of the central AGN with the galaxy's sur- 
rounding hot gaseous atmosphere. Chandra's high 
angular resolution permits the identification of 
outflow shocks, cavities, and buoyant bubbles that 
directly trace the flow of matter, energy and mo- 
mentum away from the nucleus and into the sur- 
rounding environment (e.g. see Young et al. 2002, 
Forman et al. 2005 for M87 in the Virgo Cluster; 
Fabian et al. 2003 for NGC 1275 in the Perseus 
Cluster; McNamara et al. 2000, Nulsen et al. 2005 
for Hydra A). 

It is now widely recognized that understanding 
the nature of these AGN outflows is fundamen- 
tal to further advances in the understanding of 
the evolution of structure (galaxies, groups and 
clusters) and their hot gas environments, the in- 
terstellar (ISM) and intracluster (ICM) medium. 
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Energy transfer to the ICM through episodic nu- 
clear outbursts from a central supermassive black 
hole (Binney & Tabor 1995; Owen et al. 2000) 
is likely necessary to resolve the order of mag- 
nitude discrepancy between the small amount of 
cool gas measured in the cores of rich clusters (e.g 
Tamura et al. 2001; Peterson et al. 2003 and ref- 
erences therein) and that predicted by standard 
cooling-flow models (Fabian 1994). Simulations 
demonstrated that the physical mechanisms re- 
sponsible for the transfer of energy from the re- 
sultant relativistic plasma jets, shocks, buoyant 
bubbles, and sound waves to the surrounding gas 
are complex and strongly dependent on the age of 
the outburst (e.g. Heinz et al. 1998, Churazov 
et al. 2001; Quihs et al. 2001; Briiggen & Kaiser 
2002; Reynolds et al. 2001, 2005; Ruszkowski 
et al. 2004; Fabian et al. 2005; Heinz & Chu- 
razov 2005; Omma et al. 2004). Young flows are 
characterized by shocks driven into the ambient 
medium by the outburst that heat the gas locally 
(e.g. Heinz et al. 1998; Omma et al. 2004 and ref- 



erences therein). After the outflow matures, fluc- 
tuations in the gas ( "ripples" ) passing over buoy- 
ant bubbles may induce Richtmyer-Meshkov insta- 
bilities that might dominate the transfer of bubble 
enthalpy to the ambient gas (Heinz & Churazov 
2005). 

Recent studies indicate that feedback from a 
central AGN may be fundamental to our under- 
standing of galaxy evolution in general. Sim- 
ulations, using the ACDM cosmology, suggest 
that nongravitational heating of the ISM by out- 
bursts from the central, supermassive black hole 
in galaxies is needed to correctly reproduce the 
observed shape of the galaxy luminosity function 
and massive galaxy colors at low redshift (Granato 
et al. 2004; Croton et al. 2006 and references 
therein). Best et al. (2006) used FIRST and 
NVSS data to measure the fraction of galaxies 
that are radio-active above a given 1.4 GHz radio 
luminosity as a function of that luminosity and 
central black hole mass. They interpreted their 
result as a probablistic estimate of the time a 
galaxy, with given black hole mass, is active above 
a given radio brightness, and combined this result 
with the Birzan et al. (2004) correlation between 
the current central radio luminosity and the me- 
chanical luminosity of outbursts to estimate the 
total energy available from AGN to heat the am- 
bient ISM. Best et al. (2006) argued that these 
data showed that heating from AGN could bal- 
ance cooling in typical elliptical galaxies, and that 
most of the AGN heating in these galaxies occurs 
when the radio luminosity of the central AGN is 
low, Li.4GHz < lO'^^ergss^^ Hz^^. Thus outflows 
may be important in systems with low, as well as 
high, luminosity AGN. 

In this work we use a 54.4 ks Chandra obser- 
vation to study nuclear outburst activity in one 
such galaxy, NGC 4552 {a = 12^ 35"^ 39. 8", S = 
12°33'23", J2000), an elliptical galaxy located 72' 
east of M87 in subcluster A of the Virgo clus- 
ter. In a companion paper (Machacek et al. 2006, 
hereafter called Paper I) we use the same Chan- 
dra observation to measure the properties of hot 
gas in the outer regions of NGC 4552, showing 
that these properties are consistent with that ex- 
pected for gas being ram-pressure stripped due to 
NGC 4552 's supersonic motion through the Virgo 
ICM. NGC 4552 is a member of the large sample 
of nearby elliptical galaxies used to establish cor- 



relations between X-ray emission and other ISM 
tracers, so that its global X-ray properties are well 
studied (for a brief review, see Paper I and refer- 
ences therein). 

Studies of the nuclear properties of NGC 4552, 
across a wide range of wavelengths, strongly sug- 
gest the presence of a weak AGN at its center. 
From optical line ratios. Ho et al. (1997) clas- 
sify NGC 4552 as a LINER/HH transition object, 
possessing spectral properties of a low luminosity 
AGN contaminated by nearby star-forming HII re- 
gions. However, the classification was highly un- 
certain due to the weakness of the emission lines. 
HST spectra of the nuclear region by Cappellari 
et al. (1999) showed a faint broad emission com- 
ponent suggesting the presence of an AGN and 
favoring a LINER classification for NGC 4552 's 
nucleus. Since the broad emission line component 
was observed in both the allowed and forbidden 
transitions, the nucleus was considered borderline 
between LINER types 1 and 2. Short-term vari- 
ability of the nucleus is a distinguishing feature of 
AGN activity. UV variability has been reported 
for NGC 4552 on timescales of several months 
(Maoz et al. 2005) to years (Rcnzini et al. 1995; 
Cappellari et al. 1999 ). Although not preclud- 
ing a compact circumnuclear starburst, since vari- 
ability in the UV on timescales < 1 year is not 
expected from young stars, the variable fiux pro- 
vides a lower limit on the AGN contribution in 
this waveband (Maoz et al. 2005). 

Radio continuum measurements taken with the 
VLA at 8.4 GHz (Filho, Barthel & Ho 2000) and 
with the VLBA at 5 GHz (Nagar et al. 2002) show 
a ~ 100 mjy core radio source, unresolved on milli- 
arcsecond scales, at the galaxy center. Variabil- 
ity at the ~ 20% level is observed in the 15 GHz 
and 8.4 GHz radio flux on ^ year timescales (Na- 
gar et al. 2002). Properties of the radio spectrum 
indicate a nonthermal origin for this radio emis- 
sion, most likely associated with accretion onto a 
central supermassive black hole. The measured 
brightness temperature, Tb ^ 2x 10^ K at 5 GHz 
(Nagar et al. 2002) is too high for this emission to 
be due to a nuclear starburst (Tb ^ iC — 10^ K ) 
or collection of supernova remnants (Tb < 10 K ). 
Furthermore, the shape of the radio emission spec- 
trum is too flat (a ^ 0) to be associated with a 
collection of supernova remnants, whose spectral 
indices are typically -- 0.4-0.7. Filho et al. (2004) 



attribute the flatness of the observed radio spec- 
trum to ADAF-Hke accretion onto a central black 
hole combined with synchrotron emission from 
parsec scale jets or outflows. This is supported 
by the 5 GHz observation of two ~ 5 milliarcsec 
diameter, extended features to the east and west 
of the core, suggestive of parsec scale jets (Na- 
gar et al. 2002). From the correlation between 
the black hole mass and central velocity dispersion 
(Ferrarese & Merritt 2000; Gebhardt et al. 2000), 
the mass of the central black hole in NGC 4552 is 
inferred to be 4 x IO^Mq. In recent Chandra X-ray 
observations, Filho et al. (2004) identified a hard 
X-ray source, whose X-ray spectrum is well repre- 
sented by a F = 1.51 power law component plus 
a kT = 0.95 keV Raymond-Smith thermal compo- 
nent with 0.5 — 2keV and 2 — lOkeV luminosities 
of logLx = 39.2 and 39.4 ergss"^, respectively, 
that is coincident with the core radio source in 
NGC 4552. This is also consistent with the pres- 
ence of a weak AGN. 

Our discussion is organized as follows: In §2 
we briefly review the observations and our data 
reduction and processing procedures. In §3 we 
present the background subtracted, exposure cor- 
rected image of the central region of NGC 4552 
showing two bright ringlike features within 1.3 kpc 
of the galaxy center. We model the spectral prop- 
erties of these rings, in the central 1.3 kpc of the 
galaxy, and use the X-ray surface brightness pro- 
file to estimate the density jump across the rim 
of the rings. We show that these rings are con- 
sistent with shocked gas driven outward by recent 
nuclear activity. In section §4, the spectral prop- 
erties of the nuclear region are briefly discussed 
for comparison with earlier work. In §5 we sum- 
marize our results. Unless otherwise indicated, all 
errors correspond to 90% confldence levels and co- 
ordinates are J2000. Taking the distance to the 
dominant elliptical M87 as representative of the 
distance to subcluster A of the Virgo Cluster con- 
taining NGC 4552, the luminosity distance to the 
cluster is 16.1 ± 1.1 Mpc (Tonry et al. 2001) and 
1" corresponds to a distance scale of 77 pc. 

2. OBSERVATIONS AND DATA RE- 
DUCTION 

We used a 54.4 ks observation of the Virgo el- 
liptical galaxy NGC 4552 taken with Chandra on 



2001 April 22 using the Advanced CCD Imaging 
Spectrometer array (ACIS, Garmire et al. 1992; 
Bautz et al. 1998) with ACIS-S (chip S3) at the 
aimpoint. A detailed discussion of the observa- 
tion and data analysis are presented in Paper I. 
We briefly review our procedures here. The data 
were analyzed using the standard X-ray process- 
ing packages, CIAO version 3.1, FTOOLS and 
XSPEC version 11.2. Filtering removed events 
with bad grades (1, 5, and 7) and those with sig- 
nificant flux in border pixels of the 5x5 event 
islands (VFAINT mode), as well as any events 
falling on hot pixels. The data were reprocessed 
and response files created using the most recent 
gain tables and instrumental corrections, includ- 
ing corrections for the time-dependent declining 
efficiency of the ACIS detector due to the buildup 
of contaminants on the optical filter (Plucinsky 
et al. 2003) and for the slow secular drift (tgain) 
of the average PHA values for photons of fixed 
energy. ^ Periods of anomalously high background 
(flares) were identified and removed, along with 
times of anomalously low count rates at the be- 
ginning and end of the observation, resulting in 
a useful exposure time of 51,856 s. Point sources 
were identified in the S3 chip in the 0.3 — lOkeV 
energy band using a multiscale wavelet decompo- 
sition algorithm set with a 5a detection thresh- 
old. The resulting 132 source identifications were 
excluded from the spectral and surface brightness 
analyses of NGC 4552 here and in Paper I. 

3. NUCLEAR OUTFLOWS IN NGC 4552 

In Figure 1 we show the 0.7—1 keV diffuse emis- 
sion from the two ring-like features in the cen- 
tral region of NGC 4552. Point sources were re- 
moved from the image and the point source re- 
gions filled with the local average emission level 
using CIAO tool drnfilth. In the right panel of Fig- 
ure 1 we subtract the 0.7—1 keV image smoothed 
with a 5". 4 Gaussian kernel from the same image 
smoothed with a 1" Gaussian kernel to produce 
an unsharp masked image of the same region, in 
which these rings are clearly visible. The rings 
are approximately circular and of equal size with 
radii ^ 0.85 kpc (11") and farthest projected dis- 
tance from the nucleus of NGC 4552 of ^ 1.3 kpc 



see Vikhlinin et al. in 
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Fig. 1. — {(I'Sft) Chandra image of the 0.7 — 1 keV diffuse emission from the inner rings in the central region of NGC 4552. 
1 pixel = 0."49 X 0".49 (right) Unsharp masked image of the same region, where we subtract the 0.7 — IkeV diffuse image 
smoothed with a 5". 4 Gaussian kernel from the same image smoothed with a 1" Gaussian kernel to highlight the rings. 



(~ 17"), creating in projection the 'hourglass' ap- 
pearance noted by Filho et ah (2004). Such cav- 
ities have been shown in other eUiptical galax- 
ies, such as M87 (Young et al. 2002; Forman 
et al. 2005), to chronicle past episodes of nuclear 
activity from a central AGN. 

3.1. Gas Temperature in the Rings 

We fit the spectrum of the rings in an an- 
nular region R, shown in Figure 2 and listed in 
Table 1, with inner and outer radii of 11" and 
17", respectively. The local background region 
for the rings (also listed in Table 1) was chosen 
concentric to and immediately outside the source 
to remove contamination from galaxy emission at 
larger radii as well as the Virgo ICM. We use 
an absorbed APEC thermal plasma model (Smith 
et al. 2001) corrected for absorption using Wiscon- 
sin photo-electric cross-sections (Morrison & Mc- 
Cammon 1983), with the absorbing column fixed 
at the Galactic value of 2.59 x 10^° cm~^ (Dickey & 
Lockman 1990), to fit the resulting 1840 net source 
counts over the 0.3 — 3keV energy range, where 
the source count rate was significantly above back- 
ground. We allowed the hydrogen column to vary 
and found no suggestion for increased absorption 
above Galactic. We checked that there were no 
significant differences between the spectra for the 
northern and southern rings by fitting the north- 
ern and southern halves of region R separately. 




Fig. 2. — Background subtracted Chandra image of the 
0.5 — 2 keV diffuse emission from the inner rings in the cen- 
tral region of NGC 4552 with spectral region R superposed. 
1 pixel = 0. "49 X 0".49. 



We found the temperature and abundance of the 
northern (southern) ring to be 0.64 ± 0.03 keV 
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Zq (> 0.5 Zq), respec- 



tively, for fixed Galactic absorption. Since the 
spectral fits for the northern and southern rings 
agreed within errors, we combined these regions 
into the single annular region R to improve statis- 
tics. 



The best fit APEC model spectral parameters 
for gas in tlie rings from this combined annular 
region (Region R) are listed in the first row of Ta- 
ble 2. We find the temperature and abundance 
of gas in the rings to be kT = 0.64 ± 0.02 keV 
and A = OJtoiZQ (x^/doi = 46/55), respec- 
tively. Thus the temperature of gas in the rings 
(region R) is a factor 1.4 higher than that found for 
the surrounding (0.431q'q2 kcV) galaxy gas (Pa- 
per I). This is in contrast to the properties of a 
cold front, where the brighter region is cooler than 
its less bright surroundings (Viklinin et al. 2001; 
see also Paper I and references therein) or to the 
cool, bright rims of highly evolved remnant cavi- 
ties from nuclear activity, such as those found in 
Perseus (Fabian et al. 2003) and in M87 (For- 
man et al. 2005), but is similar to the prop- 
erties observed in the X-ray edges produced by 
shocks driven into the ambient medium from re- 
cent AGN outbursts, such as in Hydra A (Nulsen 
et al. 2005a, 2005b). 

We verified that the temperature rise in the 
rings cannot be due to unresolved X-ray bina- 
ries in our data by adding a 5 keV bremsstrahlung 
component to our spectral model, as in Paper I, 
and refitting the spectrum over the 0.3 — 5 keV 
energy range. The resulting best fit temperature 
(0.63^gQ2 keV) for the thermal component in this 
two component model agrees with the single com- 
ponent APEC model results to within 2%. The 
abundance of the gas {A = 0.91q'3 Zq) is also in 
agreement, within the 90% CL model uncertain- 
ties, with that found using a single APEC model. 
The 0.5 — 2keV X-ray luminosity of region R is 
dominated by emission from galaxy gas, with the 
bremsstrahlung (unresolved X-ray binary) compo- 
nent contributing < 5% to the 0.5 — 2keV lumi- 
nosity in that region. This is in agreement with 
expectations for the residual contribution of unre- 
solved X-ray binaries obtained from a direct inte- 
gration of the average XLF for LMXB's (Gilfanov, 
2004) below our 0.5 — 8keV point source detection 
threshold (4.8 x 10^^'"' crgss^^ cm^^ ) for region R. 
Thus the effects of unresolved X-ray binaries on 
the measurement of the properties of the diffuse 
gas in the rings that arc important for our analy- 
sis are small. 

Although the uncertainties are large, the mctal- 
licity in the rings is consistent with that found 
in Paper I for the outer regions of the galaxy. 
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Fig. 3. — The 0.5 — 2keV surface brightness profiles across 
the northern (open squares) and southern (filled circles) 
rings taken in angular sectors, whose vertices are positioned 
at the nucleus and angles extend from 64° to 121° for the 
northern rim and from 257° to 349° for the southern rim 
(measured counterclockwise from the west in Figure 1) 



Thus we find no evidence for a strong metallic- 
ity gradient across the rings. For completeness we 
list in Table 2 the metallicity dependence of the 
0.5 — 2keV X-ray emissivity for gas in the rings 
over the abundance range (0.4 — 1.0 Zq) of inter- 
est for NGC 4552. 

3.2. Gas Density in the Rings 

As shown in Figure 3, the 0.5 — 2keV surface 
brightness profiles across the northern and south- 
ern rims, taken in an angular sector centered on 
the nucleus extending from 64° to 121° for the 
northern rim and 257° to 349° for the southern 
rim (measured counterclockwise from west in Fig- 
ure 1), agree within errors. Thus we use the sur- 
face brightness profile across the northern rim as 
representative of the ring system. We use the 
multivariate minimization technique (Markevitch 
et al. 2000; Vikhhnin et al. 2001) to determine 
the electron density inside the rim of the northern 
ring from the surface brightness profile shown in 
Figure 4. 

We assume a spherically symmetric power law 
distribution for the electron density for r < ri of 
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Fig. 4. — (left) The 0.5 — 2keV surface brightness profile as a function of distance from the center of NGC 4552 toward the 
north across the galaxy's leading edge from Paper I. The vertical dashed lines denote ri = 1.3 kpc and r2 = S.lkpc, the position 
of the outer radius (rim) of the northern ring and the outer leading edge of the cold front, respectively. The solid line denotes 
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the model results for the surface brightness for power law density models rig oc r " within the galaxy with a = 0.1 

r < 1.3 kpc (inside the rim of the ring) and a = 1.3_o 2 ^°^ ^■•^ < r < 3.1 kpc (between the rim of the ring and the leading edge 

of the cold front), (right) The electron density as a function of radius from the center of NGC 4552 for the above model. 



the form 



n 



(1) 



where ri is the outer radius(rini) of the ring, 
and ririm and ai are the normalization and power 
law index, respectively, for the electron density of 
NGC 4552 inside the rim of the northern ring. The 
density model for NGC 4552's gas outside the ring 
ri < r < r2, determined by our analysis of the 
northern leading edge in Paper I, is given by 



«c2 



^2 



(2) 



where r2 = 3.1 kpc, 02 = 1.3 ± 0.2, and 71^2 = 
0.01 cm^^ . Since NGC 4552 is much smaller than 
the projected 71'. 6 separation between M87 and 
NGC 4552, the contribution of the cluster ICM to 
the surface brightness for radii r < 3.1 kpc (in- 
side the galaxy) is nearly constant, varying by 
< 2%. We subtract this ICM contribution from 
the surface brightness profile within NGC 4552 
(r < 3.1 kpc) and fit the remaining galaxy emis- 
sion, taking the position Trim of the rim of the ring, 
the galaxy electron density power law indices ai, 



a2, and the discontinuity across the rim as free 
parameters. The discontinuity across the rim is 
given by 

^rim^i-im 



J 



A„ 



(3) 



lout "-out • 

where ririm and Arim denote the electron density 
and emissivity inside the rim of the northern ring, 
and riout and Aout denote the corresponding prop- 
erties of galaxy gas outside the rim. 

We find a best fit position for the disconti- 
nuity across the rim shown in Figure 4 of ri = 
1.30 kpc, coincident with the rim (outer radius) 
of the northern ring in Figure 1, and slope ai — 
0.1 ± 0.3 for the electron density at r < ri. The 
slope of the electron density distribution for radii 
r > ri (outside the rings) was a2 = 1-3 ± 0.2, 
consistent with our previous fit in Paper I for the 
galaxy gas at large radii. Assuming no strong 
abundance gradients in NGC 4552 and correcting 
for the 20% difference in the galaxy gas cooling 
functions, caused by the difference in gas temper- 
ature between gas at large radii (0.43 keV, Paper 
I) and gas in the rings (0.64 keV for region R), we 
infer the density ratio across the rim of the north- 



ernringtobenrim/^out = l-4to'2- We extrapolate 
eq. 2 to r = 1.3 kpc, the location of the rim of the 
northern ring, to evaluate riout at the discontinu- 
ity. Then using the density discontinuity at the 
rim, we find the electron density just inside the 
rim to be ririm = 0.05 cm~'^ . In the right panel 
of Figure 4, we plot the resulting normalized two 
power law density model for gas in NGC 4552 and 
the /3-model for the surrounding Virgo ICM as a 
function of distance from the center of NGC 4552. 

3.3. Shocks from Recent Nuclear Activity 

Since the temperature (fcTrim = 0.64 kcV) and 
density (ririm = 0.05 cm~^ ) of gas in the in- 
ner rings are both greater than that for galaxy 
gas outside the rings (fcTout ^ 0.43 keV, riout ~ 
0.03cm~^) at ri =: 1.3kpc, the properties of the 
rings are qualitatively consistent with resulting 
from a shock. For a monatomic ideal gas with 
adiabatic index 7, we can use the density dis- 
continuity Prim/Pout — "-rim/"-out in thc Rankiuc- 
Hugoniot shock conditions (Landau & Lifschitz 
1959) across the shock to estimate the shock speed 
(Mach number Mi) 



(7+l)M2 



nout (7 - ^)Mi + 2 
and predict discontinuities in temperature 
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and pressure 
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We can obtain an estimate for the density jump 
^rim/"-out from thc surface brightness discontinu- 
ity at thc position of the rings. However, since 
the shock front is narrow, the averaged, mea- 
sured surface brightness discontinuity will under- 
estimate the actual density discontinuity at the 
narrow shock front, and thus also the inferred 
shock strength and temperature rise. From our 
fit to the inner surface brightness edge in §3.2, 
we infer a lower limit on thc density discontinuity 
of 1.4^q'2. From Equation 4 this implies a shock 
Mach number Mi > 1.3. For this Mach number, 
the temperature and pressure increases predicted 
by Equations 5 and 6 are Trim /Tout ^ 1-3 and 



Prim/Pout ^ 1-8, in excellent agreement with our 
measured ratios of 1.48J^q}3 and 2.1^q'5, respec- 
tively. Thus thc bright rings in Figure 1 are shocks 
at thc interface between the surrounding galaxy 
gas and two expanding cavities. Numerical simu- 
lations have predicted that such features would be 
formed by high-momentum, subrelativistic bipolar 
nuclear outflows (see e.g. Reynolds et al. 2001; 
Omma et al. 2004 and references therein). 

An estimate for the age of these cavities is 
given by tage ~ R/cs where R is the mean ra- 
dius of the cavity and Cs is the sound speed in 
the ambient medium (Birzan et al. 2004). For 
R ~ 0.85 kpc and Cs ~ 340 km s^^ , the speed of 
sound in 0.43 keV galaxy gas, we find the cavities 
to be ^ 2.4 Myr old. A lower limit on the total 
mechanical energy -Emoch carried in the outflow is 
given by the work 2pgV needed to evacuate the 
two observed cavities, where pg is the mean pres- 
sure of the unshocked galaxy gas at each cavity's 
center and V is the volume of the cavity. Assum- 
ing spherical symmetry for each expanding cav- 
ity and extrapolating the density model in Equa- 
tion 2 for unshocked galaxy gas outside the rings 
to the center of each cavity, we find the sum of 
the mechanical energy for the two cavities to be 
£'moch > 1.2 X 10^^ ergs. This gives a lower limit 
on the average mechanical power for the outburst 
of Tmoch > 1-5 X 10**^ ergss^^ . 

We use a simple, spherically symmetric, point 
explosion shock model to obtain a second estimate 
for thc properties of the outflow. The details of 
this model and its uncertainties are discussed in 
Nulsen et al. (2005a, 2005b). We briefly sum- 
marize that discussion here. The model shock is 
produced by a sudden deposition of energy in the 
central cell of the simulation. The model then uses 
gas hydrodynamical equations, including gravity, 
to numerically propagate the shock radially out- 
ward through the galaxy gas, assuming thc un- 
shocked galaxy atmosphere is isothermal, spheri- 
cally symmetric with a power law density distri- 
bution, and in hydrostatic equilibrium. Assum- 
ing the gas in NGC 4552 is in hydrostatic equi- 
librium, we use our fits to the outer region (OA 
in Paper I) of NGC 4552 as representative of the 
temperature (fcT = 0.43 keV) and density distri- 
bution (oc r~^'^) of the initially unshocked galaxy 
gas. We then compute surface brightness profiles 
from the model, using the Chandra response in 



NGC 4552 data vs models 













^+ 


T ^ T-^T T— 


- 


n 


y—^ 


^'^-^-ftF^*^ 


- 


K 




V 


L-, 


,. u 


7 




^'"\ 


- 






^^^ : 


0) 






^^^^ - 


X 






^^^ 


a 






^^^ 


S^ 








CO o 


7 






C/1 ^ 






- 








* 








>_ 



r (kpc) 



> o 



0.5 



1.5 



r (kpc) 



Fig. 5. — ('e/i) Spherical shock model fits to the surface brightness profile across the inner edge for a preshocked isothermal 
(0.43 keV) r~^'^ density distribution inside the galaxy. Upper, middle and lower lines represent shock Mach numbers 1.9, 1.7 
and 1.6, respectively, (right) Temperature profile (solid line) for the Mach 1.7 model shock shown in the left panel. Dashed line 
denotes the emission-measure-weighted temperature profile for lines of sight through the galaxy. 



the 0.5 — 2 keV energy band, determined using 
XSPEC and an absorbed APEC model for tem- 
peratures determined from the shock model, as- 
suming a preshock temperature of kT = 0.43 keV 
and detector responses appropriate to this obser- 
vation. Since the shock weakens as it hydrody- 
namically evolves, and the initial conditions are 
self-similar, we scale the model flow radially to the 
position of the observed shock. We then match 
the observed surface brightness profile to that of 
the model to determine the shock parameters. Al- 
though this point explosion model will not repro- 
duce the detailed (bipolar) geometry of the out- 
flows, the model does provide a good descrip- 
tion for the shape of the surface brightness profile 
across the rim and thus a reasonable measure of 
the Mach number, and thus temperature jump, of 
the shock (Nulsen et al. 2005b). 

In the left panel of Figure 5 we see that the sur- 
face brightness profile across the rim of the north- 
ern ring in NGC 4552 is well represented by a 
shock model with Mach number Mi — I.TJ^q'^ and 
shock radius of 1.29 kpc. The predicted, depro- 
jected temperature profile across the shock front 
for the Ml = 1.7 model is shown as the solid line 
in the right panel of Figure 5. The radial extent 



of region R (0.9 < r < 1.29 kpc), used to measure 
the temperature inside the rings, is much larger 
than the radial width of the shock. Our mea- 
sured temperature ( 0.64 ± 0.02 keV from Table 
2) represents an average over this region and is in 
good agreement with the model results. Projec- 
tion effects also tend to reduce the observed tem- 
perature, since any given line of sight through the 
shock also passes through surrounding unshocked 
galaxy gas of lower temperature. The dashed line 
in the right panel shows the predicted emission- 
measure- weighted temperature profile Tproj for our 
Mach 1.7 model, taking into account the contribu- 
tion from cooler preshocked galaxy gas along the 
line of sight. The emission- measure- weighted tem- 
perature profile rises sharply at the shock bound- 
ary, peaking at ~ 0.6 keV, and then slowly de- 
creases to ~ 0.55 keV at r = 0.9 kpc, the inner ra- 
dius of region R. This predicted profile is shifted 
up or down by ~ 0.03 keV due to uncertainties 
in the temperature of the preshocked galaxy gas. 
Since the model Tproj includes the contribution 
of unshocked galaxy gas at radii larger than that 
of the rings, that we subtracted as a local back- 
ground in our previous fits to region R, we refit 
the spectrum for the rings in region R, for com- 



parison with the model, using a 6'. 7 x 2'. 5 rect- 
angular region (region VN in Paper I) centered 
at (a = 12'*35"37.7^ S = 12°36'16.5") with ori- 
entation angle of 319°, comprising approximately 
25% of the area of the S3 chip, as the local back- 
ground, to subtract only the Virgo ICM and par- 
ticle backgrounds from the data. We find a mean 
temperature /cTproj = 0.61 ± 0.02 keV and abun- 
dance ^ = 0.6 ± 0.15 (xVdof = 51.8/55) for fixed 
Galactic absorption. This is in excellent agree- 
ment with the shock model prediction (Figure 5, 
just behind the shock). We verified that this result 
was not biased by our choice of local background 
by refitting the spectrum using blank sky back- 
grounds appropriate for the date of observation 
and instrument configuration.^ When we used a 
second APEC component to model the contribu- 
tion of the Virgo cluster gas, with temperature 
2.2 keV and abundance 0.1 Zq fixed at the "best- 
fit" Virgo Cluster values from Paper I, but with 
the normalization free to vary, we found no signif- 
icant differences. 

The bipolar geometry of the outburst cavities, 
seen most clearly in the right panel of Figure 
1, makes it difficult to probe the temperature 
profile inside the rims. To investigate this re- 
gion, we extract a combined spectrum from el- 
liptical regions in the northern (southern) ring 
centered at a = 12''35"39.87^(5 = -hl2°33'31.7" 
(a = 12^35™39.53",(5 = +12°33'15.97") with 
semi-major, minor axes and position angle (a, 
b, pa) of a = 0.47kpc,& = 0.27kpc,pa = 9.5° 
(a = 0.52kpc,6 = 0.20kpc,pa = 21.3°), respec- 
tively. We use these elliptical regions rather than 
a simple annular region inside the rings in order to 
avoid the region just outside the nucleus where, in 
projection, the bipolar ring rims appear to inter- 
twine. We model the spectrum in XSPEC with an 
absorbed two component APEC model, fixing the 
temperature and abundance of one APEC com- 
ponent at the the temperature and abundance of 
the Virgo Cluster gas, as above, and hydrogen ab- 
sorption at the Galactic value. We find a mean 
temperature and abundance for galaxy gas aver- 
aged along the lines of sight passing through the 
interior of the cavities of fcTproj = 0.65lQ'o4keV 
and A = O.6t^;2^0 (xVdof = 43/34), respec- 



^see http://cxc.harvard.edu/contrib/maxim/acisbg and Pa- 
per I 



lively. This agrees within 90% confidence with the 
mean temperature and abundance (A:Tpi-oj,rim = 
0.61 ± 0.02 keV, A = 0.6 ± O.I5Z0) measured in 
the annular ring region R. 

Thus we do not observe the < 10% decline in 
temperature behind the shock that might have 
been expected from the model temperature profile 
in the right panel of Figure 5. This is not surpris- 
ing for several reasons. First, both the ring and 
cavity spectral regions average over a broad range 
in radius. Region R, with inner and outer radii 
of 0.9 and 1.3 kpc, averages over 30% of the radial 
distance from the shock rim to the nucleus, thus al- 
ready averaging over the region where the temper- 
ature decline due to adiabatic expansion of gas be- 
hind the shock front should be greatest. Similarly 
the elliptical regions for the cavities also sample a 
broad range in radius, that is further complicated 
by projection effects due to the bipolar geometry 
of the outflows. Secondly, the model evolution of 
the gas temperature behind the shock front is sen- 
sitive to the density and temperature profile of the 
unshocked atmosphere. The unshocked gas den- 
sity in galaxies typically follows a /3-model profile, 
flattening in the core. Thus it is strictly a power 
law only locally near the rim. Also galaxy gas near 
the nucleus of NGC 4552 may not be isothermal. 
Finally, the outflows are likely still driven by the 
ram-pressure of jet plasma from the initial AGN 
outburst, that has effectively evacuated X-ray gas 
from a 'cocoon' inside the shocked rims, rather 
than by a point explosion model. This would also 
affect the observed thermal properties of gas av- 
eraged over lines of sight passing through the 'co- 
coon'. In particular, if the cavities are nearly de- 
void of X-ray gas, the emission will be dominated 
by the front and back edges of the cavities. Addi- 
tional discussion of these model uncertainties can 
be found in Nulsen et al. (2005b). 

Nonetheless, the model does fit the surface 
brightness profile in the vicinity of the shock well 
(see Fig. 5), such that the shock Mach number 
is largely insensitive to the details of the model. 
Consequently, the point explosion model estimates 
for the total mechanical energy and mean me- 
chanical power carried in such outbursts are un- 
likely to err from their true values by more than 
a factor ^ 2 (Nulsen et al. 2005b). Adopting 
a factor two as a conservative estimate for these 
model uncertainties, we now use the shock model 



fit to our data to determine the total mechani- 
cal energy and mean mechanical power required 
for the nuclear outburst observed in NGC 4552. 
The elapsed time since the nuclear outburst (shock 
age) in the Mach 1.7 spherical shock model is 
1.3 Myr. Since the simplified spherical model ge- 
ometry overestimates the volume swept out by the 
observed shocks, the model predictions for the to- 
tal mechanical energy i^moch carried in the shocks 
and thus also the mean mechanical power Lshock 
of the outburst are upper bounds on the proper- 
ties of the observed outburst in NGC 4552. We 
correct for this by including a geometrical fac- 
tor /v < 1, where /v is approximately the ra- 
tio of the volume filled by the observed cavities 
to that of the spherical model. For a Mach 1.7 
shock, we find i?mech = 2.4 x 10^^/v ergs and 
ishock — 5.8 X lO^'^/vergss^^ . Assuming the 
observed cavities are well approximated by two 
spheres of radius R = 0.85 kpc, we find the frac- 
tion of the model volume filled by these caviti- 
ties is /v ^ 0.6, such that the total mechan- 
ical energy and mean mechanical luminosity of 
the outburst are -Emoch ^ 1-4 x 10^^ ergs and 
ishock ^ 3.3 X lO"*^ ergss"^ , respectively, consis- 
tent with our previous estimate. 

We obtain an estimate of the total X-ray cooling 
rate for the galaxy ISM by fitting a mean spectrum 
for NGC 4552 extracted from a 50" (3.8 kpc) circu- 
lar region centered on the galaxy's nucleus, using 
a two component absorbed APEC model with ab- 
sorption fixed at Galactic and one set of parame- 
ters fixed at the best fit Virgo Cluster parameters. 
We find a mean temperature (0.57 ± 0.01 keV), 
abundance (0.5 ± 0.1 Zq), and 0.2 - 2keV flux 
of i^(0.2 - 2) = 7 X 10-13 gj.gss-icm-2 for 
NGC 4552's ISM component, in excellent agree- 
ment with results {kT = 0.55^0:^8 keV, F(0.2 - 
2) = 6.8ti:J X 10-"ergss-icm-2) by Mat- 
sushita (2001) using ROSAT data. Using the 
0.1 — 14kcV luminosity as an estimate for the 
bolometric luminosity Lboi of the hot gas, we find 



^bol 



2.5 X lO^^ergss- 



If we use a more re- 



alistic two component temperature model for the 
ISM, consistent with our spectral results for the in- 
dividual regions, this estimate of the total fiux de- 
creases by < 2%. Therefore the energy carried in 
the observed nuclear outfiows (_E ~ 1.4x 10^^ ergs) 
is sufficient to balance cooling of the galaxy ISM 
for « 18 Myr. 



Energetic outbursts could either be powered by 
supernovae from a compact circumnuclear region 
of star formation, i.e. in the central ^ 100 pc (1".3 
radius) of the galaxy, or by an AGN at the galaxy's 
center. Assuming 10^^ ergs kinetic energy per su- 
pernova and ~ 1 supernova per 100 Mq of star for- 
mation (McNamara et al. 2004; Silk et al. 1986), 
we would need > 14, 000 supernova to provide 
the outburst energy of 1.4 x 10^^ ergs (for /y — 
0.6). This would require a mean star formation 
rate, over the 1.3 Myr age of the outburst, of > 
1 Mq yr"-'^ within 100 pc of the nucleus. If we make 
the extreme assumption that all star formation in 
NGC 4552 occurs in this circumnuclear region, we 
can use the total UV emission of the galaxy at 
2500A as a direct tracer of light from unabsorbed 
star formation, and the TIR emission ( IR emission 
integrated over the 8—1000 ^m band) as a tracer of 
star formation in obscured dusty regions, to place 
a firm upper bound on the energy from star for- 
mation, that is available to power the observed 
outbursts. Using NGC 4552's UV luminosity, cal- 
culated from photometric data at 2500A(NED) 
in Equation 1 from Kennicut (1998), we find an 
unobscured star formation rate of 0.09 MQyr""^ 
(0.14 M0yr--'^) for a continuous (starburst) star 
formation model. This is an order of magnitude 
below that needed to power the observed outburst. 
Note that the unobscured star formation rate in 
the inner 20 pc nuclear region, calculated using 
the UV flux measured by Maoz etal (2005) is very 
smaU, < 0.0002 - 0.0003 Mgyr"!. This is consis- 
tent with the high brightness temperature and flat 
spectral shape of NGC4552's radio emission that 
also argue against a dominant compact nuclear 
starburst or collection of supernovae remnants in 
the galaxy's core (Nagar et al 2002). 

The TIR luminosity provides a tracer of star 
formation in dusty regions, where the UV emission 
from young stars has been reprocessed to longer 
wavelengths by the surrounding dust. Since dust 
in elliptical galaxies tends to be centrally concen- 
trated, and dust has been identifled in the inner 
60 pc of NGC 4552 with HST observations (van 
Dokkum & Franx, 1995), this is likely a reason- 
able estimate for obscured star formation in the 
inner hundred parsecs of interest. We use the 60 
and 100 (Um FIR flux densities for NGC 4552 from 
NED to calculate the FIR luminosity (Wise et al. 
1993). We convert the FIR luminosity into the 
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TIR luminosity using a correction factor ( < 2 for 
F{60^J,m)/F{100^nJi) = 0.3; Helou et al. 1988). 
Then, using Equation 4 in Kcnnicutt (1998), we 
find the star formation rate for obscured regions 
in NGC 4552 to be < 0.03 Mgyr"!. Summing 
the UV and TIR star formation rates, we find 
a firm upper hmit for circumnuclear star forma- 
tion of < 0.12 — O.17M0yr~^, more than a factor 
of five below that needed to power the outburst. 
These low upper limits on current star formation 
in NGC 4552 are consistent with the 9.6 Gyr av- 
erage age of the galaxy's stellar population de- 
termined from optical line ratios by Terlevich & 
Forbes (2002). Thus the outburst is most likely 
powered by NGC 4552's central AGN. 

VLBA measurements (Nagar et al. 2002) show 
that NGC 4552 has a core radio source with two- 
sided east-west extended emission on ~ 5 milli- 
arcsecond scales. The radio spectrum is flat, with 
spectral index ^ (Filho et al. 2004), and ra- 
dio flux densities of 102.1 and 99.5 mJy at 8.4 
and 5 GHz, respectively. Integrating this radio 
flux over the 10 to 5000 MHz frequency band, 
we find a total central radio power of Lradio = 
1.55 X 10^8 ergs s~^ Filho et al. (2004) argue 
that the properties of NGC 4552's radio spectrum 
are most consistent with emission from accretion 
onto a ^ 4 X 1O^M0 black hole plus a self-absorbed 
compact jet or outflow. Our observation of shocks 
produced by outflows from a previous episode of 
nuclear activity supports this model. 

Furthermore, the mechanical luminosity, ighock '" 
3.3 X 10^^ ergs s^^ (with factor ^ 2 uncertainty), 
found from our model fit to the observed shock 
in NGC 4552, is consistent with mechanical lu- 
minosities inferred for cavities produced by AGN 
outbursts in other systems. Using 18 systems, 
ranging in size from galaxy clusters to the ellipti- 
cal galaxy M84 and containing well defined X-ray 
surface brightness depressions (cavities), Birzan 
et al. (2004) found a trend between the mechanical 
power Lmoch needed to create the X-ray cavities 
and the current central radio power Lradio of the 
source, that they parameterized as a power law 
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25±3 



1,-^radic 



^0.44±0.06 



(7) 



albeit with large scatter. The error on the me- 
chanical luminosity in Equation 7 reflects the sev- 
eral orders of magnitude uncertainty due to scat- 
ter in their data. From Equation 7 , we ex- 



pect the central radio source in NGC 4552 to 
produce outburst cavities with mechanical power 
LogLnicch = 41.8 ±3 ergs s^^ , where again the er- 
rors reflect the uncertainty due to scatter in the 
data around the trend line. This is in good agree- 
ment with the mechanical power LogLghock = 
41.5 ± 0.3 ergs s~^ found from our model fit to the 
observed shocks in NGC 4552, where the errors 
in our observed value reflect the factor two uncer- 
tainty in the shock model results. 

4. The Nuclear Region 

Finally, we look briefly at the X-ray proper- 
ties of the nuclear region by modeling the spec- 
trum from a 1".3 circular region (NS) centered at 
(a = 12'*35"39.8% S = -hl2°33'22.9") containing 
the nuclear source (see Table 1). After local back- 
ground subtraction using a concentric annular re- 
gion with inner and outer radii of 1".6 and 3". 2, re- 
spectively, the spectrum contains 1046 net source 
counts in the 0.3 — 5keV energy band. An accept- 
able fit to the data (x^/dof — 52/45) was found 
using an absorbed power law model with moderate 
hydrogen absorption {uh — l.l±0.4x 10^^ cm~^ ) 
and steep (F = 2.2 ± 0.2) photon index. The 
fits were improved and the photon index softened 
tor = 1.7±0.2 (see Table 3) by adding a sec- 
ond, thermal APEC component to the model and 
fixing the hydrogen absorption for both thermal 
and powerlaw components at the Galactic value 
{uh = 2.59 X 10^"cm~^ ). For these two compo- 
nent models, we first allow the temperature, metal 
abundance and photon index to vary, but find the 
abundance is poorly constrained. We then fix the 
metal abundance a,t A ~ 0.7 Zq from the best fit 
spectral model for the ring region R of NGC 4552 
(see the first line in Table 2), and find the fit is un- 
changed, although the errors on the photon index 
are modestly reduced. We checked the stability of 
this two component fit by allowing the absorbing 
column to vary and found nn consistent with the 
Galactic value. 

In the two component model, we find a pho- 
ton index F ^ 1.7 ± 0.1 and temperature kT = 
1.04^Q'^ykeV. The power law photon index is 
consistent with those expected for low lumi- 
nosity AGNs (Terashima & Wilson 2003; Filho 
et al. 2004). The temperature for the thermal 
component is higher than that found for gas in 
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any other region of the galaxy, suggesting that 
we may be directly observing the reheating of the 
galaxy ISM by the nuclear outburst. We find in- 
trinsic total 0.5 — 2 and 2 — lOkcV X-ray fluxes 
in this model to be 5.1 x 10^^"' ergs s^^ cm~^ and 
6.7 X 10""'^^ ergs s~^ cm^^ , corresponding to 0.5 — 
2.0 keV and 2 — lOkeV luminosities of 1.6 x 
Xo39 ergs s~^ and 2.1 x lO^^ergss"^. The power 
law component dominates the emission in both 
energy bands, contributing 75% and 97% of the 
emission in the 0.5 — 2 and 2 — 10 keV bandpasses, 
respectively. These results are in good agreement 
with the spectral results of Filho et al. (2004) 
for a circular region 2". 5 in diameter surrounding 
the nucleus with data extracted from the same 
Chandra observation (Obsid 2072) used in this 
work. 

Since the nuclear source is known to be variable 
at the ~ 20% level for the 2 and 3.6 cm radio flux 
on timescales ~ 1 yr (Nagar et al. 2002) and by 
a similar amount in the UV at 2500 and 3300A 
on timescales of ^ 2 months (Maoz et al. 2005), 
it is interesting to look for X-ray variability. We 
examined the light curve for a 2" radius circular 
region centered on the nucleus (a = 12''35'"39.8*, 
5 = -hl2°33'22.9") over the 0.3 - lOkeV energy 
range using 3ks binning, and found that, on these 
timescales, the nuclear count rate varied by < 2 cr 
(20%) during the 52 ks observation. We then 
searched for variability on any scale using an ap- 
plication of the Gregory-Loredo algorithm (1992) 
on the unbinned events time series. We only found 
possible variability in the lowest order to = 2 mo- 
ment (Log(OddsRatio) = 36.9, with the normal- 
ized odds ratio of 0.95), consistent with a slow 
(< 2a) increase in the mean signal over the du- 
ration of this observation. The Gregory-Loredo 
analysis failed to find any statistically significant 
evidence for variability on shorter timescales, i.e. 
in the to, > 2 moments. Thus we find no statisti- 
cal evidence in these data to support the claim by 
Xu et al. (2005), that was based on their "visual 
inspection of the light curve" , for X-ray variability 
on 1 hr timescales. 

5. CONCLUSIONS 

In this work we analysed a 54.4 ks Chandra 
observation of the elliptical galaxy NGC 4552 in 
the Virgo Cluster and found X-ray evidence for 



shocks in the inner region of the galaxy from 
a recent nuclear outburst. Two ring-like emis- 
sion features, consistent with bipolar nuclear out- 
flow cavities, are found in the X-ray images at 
r r^ 1.3 kpc from NGC 4552's center. The emis- 
sion measure weighted gas temperature through 
the rings (fcTproj ^ 0.61 ±0.02keV) and the shape 
of the surface brightness profile across the rim 
of the northern ring are consistent with a sim- 
ple spherical model for a Mach 1.7 shock from 
a 1.4 X 10^^ ergs nuclear outburst ^ 1 — 2Myr 
ago. The mechanical power carried in the shock, 
ishock ^ 3.3 X 10^^ ergs s~^, is consistent with 
that expected from the Lmech-^radio correlation 
(Birzan et al. 2004) between the mechanical lu- 
minosity and total radio power for X-ray sources 
showing nuclear outflow activity. One outburst of 
this magnitude every « 18 Myr could balance the 
radiative coohng of NGC 4552's hot ISM. 

The X-ray spectrum for the 1".3 circular region 
containing the nuclear point source is well fit by a 
two component power law + APEC model, with 
power law photon index F ~ 1.7 ± 0.1 and gas 
temperature kT = 1.04_q'j^-., in agreement with 
previous work by Filho et al. (2004). The higher 
temperature found for gas in the nuclear region 
compared to that elsewhere in the galaxy suggests 
that we may be may be directly observing the re- 
heating of the galaxy ISM in this region by the 
nuclear outbursts. 
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Table 1 
Spectral Analysis Regions for the Central 1.3 kpc of NGC 4552 



Region 


Type 


Shape 


Center 
RA, DEC 


Dimensions 
arcsec 


R 


source 


annular 


12 35 39.9, 12 33 23.4 


11,17 


R 


background 


annular 


12 35 39.9, 12 33 23.4 


20,31 


NS 


source 


circular 


12 35 39.8, 12 33 22.9 


1.3 


NS 


background 


annular 


12 35 39.8, 12 33 22.9 


1.6,3.2 



Note. — WCS coordinates for the centers of the regions are J2000. Di- 
mensions specified are radii for circular regions, and (inner, outer) radii 
for annular regions. Region R is also shown Figure 1. 



Table 2 
Spectral Models of NGC 4552's Rings 



kT 


A 


K 


A 


xVdof 


keV 


Ze 


10-5 cgs 


2Q-23gj.g(,jj^3g-l 




1 0.64 ±0.02 


0.7t°1 


5.8 


1.7 


45.8/55 


o.65t°:°^ 


0.4 


9.1 


1.1 


55/56 


0.64l0;03 


0.5 


7.6 


1.3 


48.8/56 


n 64+003 


0.6 


6.6 


1.5 


46.4/56 


0.64 ±0.02 


0.8 


5.1 


1.9 


46/56 


0.64 ±0.02 


1.0 


4.2 


2.4 


47.4/56 



Note. — Spectral model parameters for the rings from region R 
(see Table 1 and Figure 1), using an absorbed single temperature 
APEC model fit over the 0.3 — 3keV energy range. The spectrum 
contains 1840 net source counts. The hydrogen absorption column 
is fixed at the Galactic value (nn — 2.59 x lO^o cm~^ ). ^ The first 
row denotes the best fit model with temperature and abundance 
free to vary, while following rows hold the abundance fixed. Errors 
are 90% confidence limits. 
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Table 3 
Spectral Models of the Nuclear Region of NGC 4552 



102"cm-2 


kT 
kcV 


A 


a 


xVdof 


11±4 
2.59-^ 
2.59-^ 
<6.7 


1 04+°l-'' 

1 04+0-21 
^■"^-0.11 

-I n/i+0.13 
l-04_o.i8 


5+0-^ 
0.7/ 
0.7-^ 


2.2 ±0.2 

1 7+0.2 
l^-'-0.3 

1.7±0.1 
1.7 ±0.2 


51.7/45 
35.6/43 
35.7/44 
35.6/43 



Note. — The spectrum contains 1046 net source 
counts in region NS, specfied in Table 1, over the 0.3 — 
5keV band used for these fits. Superscript / denotes a 
fixed parameter. Errors are 90% confidence hmits. 
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